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Abstract 

It has been observed recently that a consistent LO BFKL gluon evolution leads to a 
steep growth of F2{x,Q^) for x —>■ almost independently of Q^. We show that current 
data from the DESY HERA collider are precise enough to finally rule out a pure BFKL 
behaviour in the accessible small x region. Several attempts have been made by other 
groups to treat the BFKL type small x resummations instead as additions to the conven- 
tional anomalous dimensions of the successful renormalization group "Altarelli-Parisi" 
equations. We demonstrate that all presently available F2 data, in particular at lower 
values of Q^, can not be described using the presently known NLO (two- loop consistent) 
small X resummations. Finally we comment on the common reason for the failure of these 
BFKL inspired methods which result, in general, in too steep x- dependencies as x —>■ 0. 



(To be published in Phys. Lett. B) 



In spite of the success [|I], ^ of the conventional renormahzation group "Altarelh- 
Parisi" equations (RGE) 0] much work has focused recently on replacing or modifying 
them by small x resummations. It should be stressed that, so far, this is not required by 
the experimental data. The rise of F2 at small x, for example, has been predicted by the 
GRV group P using the NLO RGE and thus is not a unique sign of the BFKL type 
of growth. By using the BFKL evolution methods of Askew et al. [§] consistently, we did 
even show in a previous paper that it is probable that a pure LO BFKL evolution 
is not compatible with the observed growth in the F2{x,Q'^) data. We confirm this here 
using more recent data from the DESY HERA collider [0, |^ . 

But this does not necessarily prove that the BFKL resummation of leading logarithms 
in X is inadequate. The original resummation ignores the quarks and uses a fixed strong 
coupling constant. Thus it is expected to be relevant only in the asymptotic small x limit 
and it could still be valid at very small values of x that are currently not accessible by 
experiment. There is also the possibility that the corresponding NLO resummations will 
change the LO behaviour and that the assumption of gluon dominance limits the region 
of applicability strongly. 

Keeping this in mind the "high energy" [s ^ leading to small x via x ~ Q'^/s) 



resummations done by Catani et al. |T0[ seem to offer a way out. They treat the BFKL 
type resummations as additions to the fixed order anomalous dimensions and Wilson 
coefficients, hoping to obtain an improved perturbative expansion that way. They also 
resummed the NLO small x contributions in the quark sector. Furthermore, they calcu- 
lated all resummations consistently by using the fcT^-factorization and showed that their 
methods do not spoil the all-order factorization of coUinear singularities. We note that 
by recalculating the LO resummations they have confirmed the known BFKL results. 

The obvious question is then of course if the "improved" expansion can match the 
success of the conventional RGE. First comparisons with F2 data using these methods 
were made by Ellis et al. . They stressed the importance of enforcing the fundamental 
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energy-momentum conservation which is generally broken by the resummations. Starting 
the evolution of fiat input densities from Qq = A GeV^, they can describe the F2 data 
fairly well. In this paper we will show that the agreement is spoiled when the starting 
scale is lowered to values typical for recent RGE fits, in order to take into account all 
presently available HERA data. 

First we will present an update of the analysis of F2 data using the LO BFKL evolution 
of the gluon. We give now a brief outline of the methods employed in , noting that the 
references for the formulae and methods used can also be found there. The BFKL equation 
concerns the unintegrated gluon density f{x,k'^), because the transverse momentum fc^ 
of the gluons in the corresponding ladder diagram is not strongly ordered as in the LO 
RGE case and cannot be integrated out. 

The connection between the integrated and unintegrated gluon is given by the equa- 
tions 

X9{x,Q') = I -j^fix,^), (la) 







f{x, k 



2^ 



dx g{x, 



(lb) 

Q2 = A;2 



91nQ2 

and in terms of the unintegrated gluon one can then write the BFKL equations as an 
evolution equation in x 

V(x,fc'2)-/(X,P) , fix,k') 



df{x,k^) _ 3a,(fc2) f dk'^ 

X ~ — rC 



dx TT ./ k'"^ 



\k'^-k^\ (Ak'^ + k^)^ 



(2) 



Predictions for F2 are then made by convoluting the BFKL gluon with the photon- 
gluon fusion quark box F^'^^ using the /cT^-factorization theorem 

F.{x,Q') = ['^f( -,k'A Fr\y,k-,Q% (3) 



k'^ 

with i = T,L denoting the transverse and longitudinal parts, respectively. Since we 
have completely ignored the quark sector and non-perturbative effects, we will add a 
background fitted to the data. Obviously the background should not dominate over the 



BFKL contribution and it also should not have a strong x dependence that would mask 
the BFKL behaviour. We use a soft Pomeron ansatz Cipx~^'^^, where the constant Cip 
is fitted separately for each bin. It is well known that the BFKL equation has the 
same double logarithmic limit as the LO RGE, when its (fixed) coupling constant as is 
replaced by the standard running one of the LO RGE. Since we want to connect to the 
LO RGE at large x, this replacement has to be made in (||). 

LO RGE gluon distributions are then used at xq = 10~^ as input via Eq. (|l^). We 
now have to face the difficulty that requires knowing /(x, fc^) down to A;^ = GeV^, 
but of course even the special LO RGE gluons g^^^^XjQ"^) used in |0 do not extend far 
enough down in to supply that information. Thus we are forced to make an ansatz 
for the infrared (IR) region, if we want to avoid cutoffs to the integral. The introduced 
running coupling also needs to be taken care of in the region of small fc^. Following Askew 
et al. the IR ansatz of is 

/(x,^^) = ^^^^/(a:,fc,2) for P < kl (4) 

At X = Xq we additionally shift the gluon 

/^«^(xo,fc2) — /^«^(xo,fc2 + A;D. (5) 

The strong coupling is frozen: as(fc^) — > ctsik"^ + k^). The ansatz is scrutinized closely 
in We just want to remind the reader of the salient points: old gluon distributions 
have to be used, since most of the more recent ones already create a steep growth of F2 
in the low x region without any additional BFKL boost. The dependence of F2 on the IR 
parameters is under control, except for k^ which has to be fixed by forcing the integrated 
gluon of @ to match the original gluon g^'^^ at medium to high Q^. For the Dp-type 
gluon used below we obtain k'^ = 0.95 GeV^ 0. The strong BFKL growth of F2 could be 
suppressed by increasing k^, but this would of course spoil the consistency with the 
RGE gluon. 

Parametrizations of our results were given in . Using the Do-tjpe gluon parametriza- 
tion and fitting the background to the new data below x = 10~^ from the HERA collab- 
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orations |^ at DESY as well as to data from the Fermilab E665 experiment |T^, we 
obtain the solid curve in Fig. |l|. The background contribution is shown as dashed curve. 
Note that at = 3.5 GeV^ and at = 6.5 GeV^ C/p would be negative and is set to 
zero by hand. It is obvious that in spite of having the freedom of fitting Cjp, the data 
can not be described by the BFKL curves. The steep growth of F2 ~ x~^'^ predicted by 
the BFKL equation is simply too strong, especially at medium to low Q^. Thanks to the 
precision of the new F2 data pure BFKL evolution is now definitely ruled out! 

Of course immediately the question arises if this can be ameloriated by using modified 
evolution equations that incorporate the successful RGE. That the NLO ROE are very 
successful indeed at describing the F2 data is made obvious by the dotted line in Fig. ^ 



which shows the Rl fit of MRS |]I3|. We are especially interested in elucidating the 
question if modified equations can match the success of the NLO RGE in the region 
of low to medium Q^. The fact that even below = A GeV^ data can be described 
using those equations has been pioneered by the GRV group ||5|, |^ and is now becoming 
universally accepted p|, ^ [l^, |I4|- This region is also expected to lead to difficulties with 
the resummed BFKL pole at A = 4 In 2, because as increases strongly in this region 
if it runs. 



We have tried the method employed by Forshaw et al. [|T5|. They limit their calcula- 
tions to the LO small x resummations of the gluon by replacing the standard anomalous 
dimension in the RGE by the resummed one and by ignoring the quarks. The (unsuc- 
cessful) results of trying to extend this method down to lower will be reported in a 



forthcoming detailed article |]T6[. Instead we concentrate here on the method developed 



by Ellis et al. |TT|. Following the work of Catani et al. ||T0| they add the known LO and 



NLO resummation corrections to the NLO RGE singlet anomalous dimensions and to the 
Wilson coefficients. Note that the corrections to the non-singlet anomalous dimensions are 
less singular than 1/x for x and need not be taken into account. We use the Mellin 
transformation with respect to x in order to express the results in terms of moments 
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f{n,Q')= / dxx-f{x,Q'), (6) 







where the MeUin moment is shifted by one compared to the usual convention (x" since 
this will allow a more compact notation in the following. 

To the fixed order expansion of the anomalous dimensions 

7(n,a,) = ^j^'\n) + (^Y j^'\n) + 0{a',) (7) 

Z TT Vz TT/ 

towers of small-x resummation corrections are now added. In LO, that is summing terms 
of the order (ctg/n)'^, only corrections to the anomalous dimensions and exist, 
where E = (q + q) denotes the singlet quark density. In NLO the corrections to 7se 
and to 7e£, are known, but those of the gluon have not yet been completed. The additional 
resummation corrections can thus compactly be written as 



Tres. 



res. 



7SS 7Sg 

7ffS 799 

gf7L 7^; + ^ 0(?) 0(?) 

with A; > 0, and the question marks signify the yet unknown entries. 
To obtain the BFKL gluon anomalous dimension 71, from 

n^asxilL); as = ^-^; x(7) = 2 V'(l) - -0(7) - "0(1 - 7), (9) 

TT 

we use a complex Newton method. For the NLO resummation obviously the choice of 
scheme is important. Since up to now completely resummed expressions are only available 
in the DIS factorization scheme, we use this scheme in the following. 

The DIS scheme is defined by Cf^ = 1 and C| = 0, thus all the resummation corrections 
are shifted into the anomalous dimensions. Then we can write 

INL = -^h2ijL,n)R{-fL), (10) 

/ 

h2{l,n) = h2{l) {1 + 0{n)) for n ^ 0, (11) 
^2(7) = — (,T34^sin(7r7)tan(7r7);' ^^^^ 



i?(7) 



r(i-7)x(7) 

r(l + 7) [-7X'(7)] 



Here the prime denotes the first derivative with respect to 7. The expressions can be 
directly calculated once 72, has been determined. 

It is a non-trivial check of the work of Catani et al. that by expanding the small x 
resummations in as/n one obtains the same constant respectively 1/n singular terms as 
in the (NLO) RGE in the — > limit. To avoid double-counting, these terms have to 
be subtracted before the additional resummations are added. Effectively, Eq. (|) is then 
replaced by 



7r, 



73 



2%8f7g 2ny7, 
gl7. 7. 
as 

IL 

n 



1 



13 as 
6 n 



(14) 

(15) 
(16) 



Here we can see clearly that the corrections to the splitting into gluons and quarks just dif- 
fer by a colour factor. It should also be pointed out, that the first terms in the expansions 
of the effective corrections are 

7, = 2C(3) f^V + 0( (-) 1, (17) 



n 



7g 



^ 2 
2n ^3 



71 
18 




(18) 



so that the correction in the quark sector leads in terms of powers of and also starts 
with the lower power of 1/n. Thus we expect the quark correction to play an important 
role despite its being formally subleading. The quark correction is also expected to vanish 
more slowly for x —>■ 1, because it does not fall off as steeply when n 00. 



It is of course entirely possible that the NLO order corrections in the gluon sec- 
tor change those expectations, since they depend on the fact that the second and third 
coefficients in the expansion of 7^ vanish. A further problem is that the added correc- 
tions violate the fundamental energy-momentum conservation, as can be easily checked 
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by setting n ^ 1. The calculation of higher resummation corrections will not solve 
this problem, since the LO corrections cannot be cancelled. Of course the original RGE 
manifestly obey energy-momentum conversation. This is not only a theoretical prob- 
lem, since the unconserved growth is very strong and leads to obviously false results. 
We will look at the simplest ways to restore the conservation of energy and momentum. 
The "hard" way is simply to subtract the troublesome values of the corrections via 
7res.(^) 7res.('^) " 7res.(l)- As has already been remarked in [jll| this does not much 



reduce the growth and it has turned out to be impossible to fit the experimental F2 data 
using this approach. Thus we will limit ourselves to the second "soft" method, multiplying 
the corrections by a factor that goes to zero for n —>■ 1. 

This conserving factor, being introduced in an entirely ad hoc way, should go to unity 
in the limit of very small n which corresponds to small x. We test the following factors: 
(1 — n), (1 — n)^ and (1 — 2n-|-n^). They have been studied previously in |T^ using higher 



starting scales Qq. Obviously these factors will inhibit the growth more strongly, since 
they weaken the singularities as can be seen by looking at the expansions (p!7|) and (plsl) . 
It should also be noted that multiplying the quark correction by the third factor leads 
to a term ~ n. This means that the correction times this factor cannot be transformed 
back into x-space on its own. Instead it needs a parton distribution that falls off at least 
with 1/n. While this is fulfilled for all reasonable parton distributions, it is still a hint 
that the parton distributions are not correctly factorized anymore. For this reason we do 
not examine even higher powers of n in the conserving factor. By checking this range of 
factors we also get an estimate of the influence of the missing higher order corrections, 
since they will also be suppressed in powers of 1/n versus powers of 

The RGE are only known to two loops and can be conveniently solved in n-space. We 
wish to modify the known NLO RGE n-space solution |TB|, but then it is necessary to 
check if the corrections can be added to the two-loop solution consistently. Since we only 
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know the running couphng to two loops, we can write the evolution equation as 



dq{n, 
da. 



^ g(n,Q^), (19) 



4n (47r)2 



using for the strong coupling (n/ is the number of active massless quark flavours) 
c^siQ') 1 A Inln(gVA^) . . 2 . ^8 

In it is shown that a consistent addition of the corrections can be accomplished by 
replacing 



We have tested that by using this method we can reproduce the results shown in |TT|. 



The dashed curve of Fig. ^ shows the strong effects of the added resummations 
compared to the solid curve of the conventional RGE. For both curves the unchanged 
Rl input distributions at Ql were used and the resummations were added using the 
conserving factor (1 — n). Obviously it is necessary to refit the input parton distributions, 
since the growth in the small x region is incompatible with the data. It is interesting 
to note the similarity of that growth to the LO BFKL one shown in Fig. ^ as dotted 
line. This is a first hint that fitting experimental data with the new method will lead to 
the same problems at small x. In order to test if the RGE plus corrections can fit the 
experimental data as well as the pure RGE, we would in principle need to fit to the same 
range of experiments as is done for example in [|T^. Even if the resummation corrections 
to all those processes were known, the computing time would still be prohibitive due to 
the complicated calculation of the corrections. 

Instead we use the following approach to save computing time: we take the optimal 



MRS Rl |T3| fit as the basis for our fit by using the same starting scale Qq, J^qcoinf = 4) 
and the same valence distributions for the quarks. Since the Rl fit describes the HERA 
data very well, we use it with the unmodified RGE to generate F2 points for our fit. Below 
X = 5 ■ 10~^ the errors and x-range used for the points correspond to those of the HERA 
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data. Above this value of x we force our fit to reproduce the Rl fit by setting an artificial 
error of one percent. In this large x region not only F2 but also the sea and the gluon are 
tested. This simulates the numerous experiments constraining the parton distributions 
at large x which are taken into account by the Rl fit. Since we do our calculations in 
the DIS scheme, we have to apply the corrections for the change of scheme to the MS 
partons of Rl. This is done at the starting scale Qq = 1 GeV^ and then the RGE in the 
DIS scheme are used to obtain the constraints. In fact the fitting program uses the same 
transformation, so that the new fitted parameters are valid in the MS scheme and can be 
directly compared to those of the MRS Rl fit. 

It turns out that the modified evolution equations cannot reproduce well the RGE sea 
in the large x region. The gluon and F2 are reproduced with reasonable accuracy. But the 
sea is less and less in agreement when the maximum power of n in the conserving factor 
is raised. This can be easily understood by noting that the product of the conserving 
factor with the expansion (|T8D does not fall off as quickly for n ^ 00 as the product 
with (^) does. Relaxing the strong constraints at large x is not a viable solution to 
the problem. This is so because, as we will see, the F2 data at small x are not well 
described by the modified evolution equations. Thus once the tight constraints in the 
large x regions are removed, the fit program adjusts the parton distributions for the small 
X region and basically ignores any weak constraints imposed at large x. This leads to 
strong deviations especially of the sea from the original Rl partons. Since here our goal is 
not to give an alternative description of the small x data, but rather to check if modified 
evolution equations can match the RGE in describing all data, we cannot avoid imposing 
the constraints at large x. 

The results of our fits are displayed in Fig. |^. The curve denoted as "no conservation" 
is the fit without conversation of energy and momentum, corresponding to the original 
ideas of Catani et al. [|1^. The one described as "shift" is a fit using the conserving factor 
(1 — n) but also shifting the strong coupling using as{Q'^) ctgiQ"^ + 1.5 GeV^). We 
note that the fit using the conserving factor (1 — 2n + n^) is not displayed. It leads to 
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an unacceptable fit which oscillates around the low x data strongly, especially at the low 
values. This confirms our earlier suspicions about possible problems with this factor. 
The "no conservation" curve shows that the unconserved growth is way too strong to be 
in accord with experimental data. 

It is also easy to see that increasing the maximum power of n in the conserving factor 
leads to the expected decrease of the growth in the small x region. Still it is obvious that 
the curve corresponding to the (1 — n) factor cannot describe the data except perhaps at 
the highest bin. The curve pertaining to the (1 — n)^ factor offers some improvement. 
This is expected due to the better cancellation of the leading term in the added quark 
resummation (p!8D for n < 1 which leads to a smaller deviation from the standard RGE. 
But the small x region below 12 GeV^ is still not fitted well. That this is caused by the 
strong BFKL like growth at small x is especially obvious at 1.5 GeV^. Here we can see 
that the fit has adjusted the parton distributions so that the curves at medium to low x 
are below the data to compensate the growth at very small x somewhat. 

Looking at the "dips" in the curves we can also clearly see that the onset of the growth 
is shifted towards lower x when the power of n in the conserving factor is increased. We 
can conclude from these fits that the effects of changing the conserving factor are sizable 
in the small x region and that no definite predictions on F2 at small x can be made for the 
time being using modified evolution equations. On the other hand the observed strong 
rise at small x could be taken as an indication that the low region will not be better 
described once all NLO terms of the resummations compatible with the two-loop RGE 
are knownQ. The BFKL like growth simply is too strong for low and spoils agreement 
with data (rather flat in x) whenever it becomes dominant. This is shown nicely by the 
curve labeled "shift" which is generated using a shifted value of a^. Obviously the fit to 
the data is much improved at low Q'^ since the BFKL growth is suppressed due to the 

artificially lowered value of as- 

^ However higher order terms introduced by a purely kinematic constraint on the gluon ladder suggest 
that the growth of the gluon could be reduced . 
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All our attempts to describe experimental data of F2 using methods inspired by the 
BFKL equations have failed. The direct application of the BFKL equation in the small 
X region is now finally ruled out by experimental data [|2|, ^, [12|. Adding the small x 



resummations to the RGE does not lead to a conflict with data, if the parton distributions 
are adjusted accordingly and provided the starting scale Ql is kept around 4 GeV^ |11|], 



i. e. if only a fraction of all available data is taken into account. But it is by now widely 
accepted that it is possible to describe data below this Ql scale using the NLO RGE. The 
modifled evolution equations fail this acid test. This conclusion can be drawn despite the 
apparent difficulty to obtain any definite predictions using these equations at all. 

A further indication of the reliability of our results is provided by a similar analysis 
of Ball and Forte pOf. They also found that deviations from the standard RGE are not 



supported by recent HERA data, although their method differs in some respects from 
ours. Notably they switch on the small x resummations at a fitted xq, they estimate the 
missing NLO gluon small x contributions by using energy-momentum conservation, and 
refit only the small x tails of given parton distributions. 

There is a common reason for these failures. At low the BFKL pole A = 4 In 2 
dominates over all the other poles since it moves strongly along the positive real axis in 
the n plane due to the growth of as- The validity of this argument is nicely demonstrated 
by the curve labeled "shift" in Fig. ^ Simply calculating the values for A at = 4 
GeV^ and at = 1 GeV^ suggests that any future method that includes this pole in the 
evolution is likely to have trouble to describe the moderate rise in x of F2{x, Q^) observed 
in the experimental data at lower values of Q'^. 
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Figure Captions 



Fig. |l] BFKL predictions using the methods of compared with recent data 0, [12 



The sohd curve is the LO BFKL prediction including the fitted background which 
is shown separately as dashed curve. Only data with x < 10~^ is fitted. The dotted 
curve shows the conventional RGE fit Rl of MRS JISl for comparison. 



FIG. The original MRS-Rl input |T3| evolved with the RGE (solid curve) and the 



RGE plus resummations (dashed curve). The resummations are added using the 
conserving factor {1 — n). For comparison the LO BFKL prediction is shown as 
dotted curve. Note that "Rl (res.)" is not fitted to the shown data 

Fig. ^ Modified evolution predictions, based on the RGE supplemented by the small x 
resummations according to Eqs. ([l^ and (|2T|) , compared with experiment (data as 
in Fig. |l]). The stars show the generated data used for the fit as explained in the 
text. The solid and dot-dashed curves are obtained with (l—n) as conserving factor, 
using a shifted in the latter discussed in the text, and the dashed curve 

is calculated with (1 — n)"^. The dotted curve is obtained without any conserving 



factor, which corresponds to the original results of |10 . 
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